Interacting damage models mapped onto Ising and percolation models 
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We introduce a class of damage models on regular lattices with isotropic interactions between the 
broken cells of the lattice. Quasistatic fiber bundles are an example. The interactions are assumed 
to be weak, in the sense that the stress perturbation from a broken cell is much smaller than 
the mean stress in the system. The system starts intact with a surface-energy threshold required 
to break any cell sampled from an uncorrelated quenched-disorder distribution. The evolution 
of this heterogeneous system is ruled by Griffith's principle which states that a cell breaks when 
the release in potential (elastic) energy in the system exceeds the surface-energy barrier necessary 
to break the cell. By direct integration over all possible realizations of the quenched disorder, 
we obtain the probability distribution of each damage configuration at any level of the imposed 
external deformation. We demonstrate an isomorphism between the distributions so obtained and 
standard generalized Ising models, in which the coupling constants and effective temperature in 
the Ising model are functions of the nature of the quenched-disorder distribution and the extent 
of accumulated damage. In particular, we show that damage models with global load sharing are 
isomorphic to standard percolation theory, that damage models with local load sharing rule are 
isomorphic to the standard Ising model, and draw consequences thereof for the universality class 
and behavior of the autocorrelation length of the breakdown transitions corresponding to these 
models. We also treat damage models having more general power-law interactions, and classify the 
breakdown process as a function of the power-law interaction exponent. Last, we also show that the 
probability distribution over configurations is a maximum of Shannon's entropy under some specific 
constraints related to the energetic balance of the fracture process, which firmly relates this type of 
quenched-disorder based damage model to standard statistical mechanics. 

PACS numbers: 46.50. +a,46.6.5.+g, 62.20.Mk, 64.60.Fr,02.50.-r 



I. INTRODUCTION 

The physics of breakdown processes that lead, for 
example, to stress-induced catastrophic failure of man- 
made and geological structures, remains an ongoing sub- 
ject of research. Stress-induced fracture of a homoge- 
neous material containing a geometrically simple single 
flaw has been studied since the work of Grifflth 01 and is 
now well understood. However, the breakdown of hetero- 
geneous structures, in which the local mechanical proper- 
ties are randomly distributed in space and/or time, con- 
tinues to present challenges despite the many advances 
over the last flfteen years The difficulty is in charac- 
terizing and quantifying the effects of interaction between 
the multitude of constituents. 

Most of the knowledge about these types of systems 
has been obtained from lattice network simulations. One 
of the most well-studied lattice network models is the 
80 year old Fiber Bundle Model (FBM) 0, H H that 
describes the rupture of bundles of parallel fibers. This 
model originally considered elastic fibers of identical elas- 
tic constant, breaking when their elongation exceeds in- 
dividual thresholds distributed according to a given un- 
correlated random distribution. A global load sharing 
rule (GLS) is assumed by which the load carried by 



a fiber is uniformly distributed to the surviving fibers 
when it breaks. Analytical results for this model have 
been obtained 0, 0, S that concern the aver- 

age load/ deformation properties ^i, iSi, iSiJ, the distribu- 
tion of avalanches or the relationship between such 
quenched-disorder based models and standard statistical 
mechanics |1L] . 

The original FBM model with global-load sharing has 
been generalized to allow for non-uniform load shar- 
ing rules, considered either as purely local load shar- 
ing (LLS), in which case distribution of avalanches and 
mechanical properties have been analytically studied 
|il2i. .1,1. .141. or as power laws of distance from the failed 
fiber, which have been studied numerically 

IHEl. The 

distribution of avalanche size s has been shown to be a 
power law s~^/^ close to macroscopic breakdown in GLS, 
and not to follow any power-law in LLS 0, 0| . 

Closely related to these models of fiber bundles in elas- 
tic interaction, many studies have addressed fuse net- 
works 11, 18, 19, 20, 21, 22,, ZililliBl, or networks of 
isotropic damage with interactions having the range of 
the elastostatic Green function . Such work 

has recently helped to understand the origin of a uni- 
versal geometric feature of crack surfaces; namely, their 
large-scale self-affinity. It is now established that three- 
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dimensional fracture surfaces in disordered brittle solids 
are self-afEne, with a material- independent r oughne ss ex- 
ponent C ~ 0.8 at large scales HHElEl 
a cross-over to C ~ 0.5 at smaller scales. 
In another load geometry, interfacial crack pinning be- 
tween two sintered PMMA plates have been reported to 
give rise to crack fronts with C ^ 0.6 |^. Recent 
work ji^, I3| has explained these large-scale roughness 
exponents using gradient-percolation theory {^^, rela- 
tionship to fracture problems first introduced in Refs. 
[3 EBl); together with estimates of the correlation- 
length divergence exponent obtained using finite-size 
scaling of lattice models in the approach to failure. Scal- 
ing laws were produced relating the roughness exponents 
of fracture or damage fronts, to the critical exponent de- 
scribing the divergence of an autocorrelation length in an 
appropriate damage model. Thus, one important objec- 
tive of the present work is to analytically extract such 
divergence exponents for a broad class of damage models 
and to show how such exponents depend on the model 
characteristics. 

The classification of lattice breakdown processes as 
critical-point phenomena is still subject to debate [flL liiil 
EilillilliililEliSlil. Local load sharing mod- 
els are usually understood as breaking through a process 
similar to a first-order transition while models with 
long-range elastic interactions or GLS are analyzed either 
as a critical-point transition licl l4ll li^ . or as a spin- 
odal nucleation process l49| . The issue depends on 
whether the accumulated damage has a correlation length 
that diverges as a power law of the average deformation 
in the approach to macroscopic failure. Numerical ev- 
idence in the literature suggests that the nature of the 
correlation length in the approach to failure depends on 
the specific model being analyzed, the damage interac- 
tion range, and the type of quenched-disorder distribu- 
tion considered. A major difficulty of such attempts to 
classify the breakdown process is that no analytic form 
for the distribution of the damage configurations at a 
given external load is available (at least that has a firm 
basis). Numerical simulations of the scahng behavior of 
avalanches are limited in size due to computational re- 
strictions which makes critical-point analysis difficult. 

It is proven here that quasi-static interacting damage 
models, that possess local breaking thresholds randomly 
quenched ab initio, can be mapped onto percolation, 
Ising or generalized Ising models with non-zero couphng 
constants for distant cells, depending on the range of the 
assumed load sharing rule and on the specific probabil- 
ity distribution of the breaking thresholds. The stress- 
induced emergent-damage distribution is proven here to 
be a Boltzmannian with a temperature (probabilistic en- 
ergy scale) that is an explicit analytical function of the 
applied external load. Having an analytical form for 
the damage-state distribution function allows a rigorous 
classification of the breakdown process into first-order or 



critical-point phase transitions belonging to various uni- 
versality classes. 

We have recently proposed a statistical theory for the 
localization of oriented fractures that emerge and elas- 
tically interact when the system has a shear stress ap- 
plied to it IHnl . I51I IH^ . The distribution of the emer- 
gent crack states was obtained using the postulate that 
the fracture arrival would maximize Shannon's entropy 
under constraints representing the energetic balance of 
the process. In the present work, we will not make this 
postulate but will instead prove its validity by direct in- 
tegration over the damage evolution. The present work 
also presents more general ranges of interaction. How- 
ever, unhke our previous work, the present analysis is for 
a purely scalar description of damage interaction. The 
interaction of real fractures in an anisotropic load, where 
microfractures present high aspect ratios, necessarily re- 
quires a tensorial elastic description as our earlier work 
|H3,E1E3 provided. 

Other analytical treatments of fracture processes in 
disordered brittle materials can be found in the litera- 
ture [ifll IH^. These approaches directly postulate both 
the appropriateness of classical Boltzmann distributions 
and the form of a free energy containing non-local long- 
range interactions. Such approaches have led to the anal- 
ysis of fracture processes as a spinodal nucleation. How- 
ever, the mechanical origin of such a free energy and its 
relation to the underlying presence of fractures and het- 
erogeneity in the system, as well as the physical meaning 
of the temperature parameter, remain unclear in these 
models. Fractures usually form irreversibly, which makes 
such direct applications of standard equilibrium statisti- 
cal mechanics questionable. Fracture theories as a ther- 
mally activated process have long been proposed. 
However, the pertinence of thermal activation in highly 
heterogeneous and stiff systems is doubtful, because the 
thermal energy is negligeable compared to the gaps be- 
tween the strain energies of the possible states of the 
system. 

Thermally induced fracture models generalized 
with quenched disorder in the rupture thresholds [HHES 
|H3,Ea|, have recently been confronted with experiments. 
Microacoustic emissions in heterogeneous materials prior 
to failure were experimentally observed to have a cumu- 
lative energy that follows a power law in time to failure, 
while the individual microacoustics events have energies 
distributed as a power law 1^ ■ In models repro- 

ducing these experimental results through time-delayed 
fiber bundles with GLS, it has been pointed out that the 
type of quenched-disorder distribution assumed for the 
individual failure thresholds directly infiuences the equiv- 
alent ternperature in thermally activated rupture models 
[HE EE EE l^Ji which makes this temperature signifi- 
cantly different from the usual temperature associated 
with molecular motion. 

In comparison with these latter works, the models de- 



3 



veloped here are quasistatic and do not include any time- 
delayed fracture processes, but their advantage is the con- 
sideration of an arbitrary range of interactions and of ar- 
bitrary distribution functions for the fracture thresholds 
(similar to Refs. [HEl, 

except that the present work 
is analytical rather than numerical) . The nature of the 
energy scale (called temperature) that enters the proba- 
bility distribution will also be clarified. 

The organization of the paper is as follows: in section 
im we introduce (and justify in the appendix) the gen- 
eral type of scalar damage models that are considered. 
In section IIIIL the probability of each damage configura- 
tion is obtained by integrating over all paths that lead 
to it. In section IIVI we estabHsh the relationship be- 
tween these configurational distributions and standard 
statistical mechanics, which allows the standard toolbox 
of statistical mechanics to be applied to damage models 
entirely based on quenched disorder. These analytical 
developments will then be utiHzed in section to estab- 
lish that fracture processes in such damage models are 
isomorphic to percolation for GLS l|V All , or to the Ising 
model for LLS dV Bp . This allows us to isolate some tran- 
sition points in these models, and to predict the nature 
of the correlation length in the approach to the tran- 
sition. We will also discuss in Section IV CI the case of 
damage models with arbitrary power-law decay of the 
interactions, and show how they are related to general- 
ized Ising models, which can themself be mapped onto 
standard ones via renormalization of the coupling con- 
stants. The results are summarized and discussed in a 
concluding section. 

II. DEFINITION OF THE DAMAGE MODELS 
CONSIDERED 

Our models reside on a regular lattice of dimension D 
(e.g., a square lattice in D = 2). Each cell of the lattice 
has a location x within the ensemble ft of cells making up 
the lattice, and has a state described by a local order pa- 
rameter (fix , where ipx = if the cell is intact and (fix = 1 
if it is broken. There is a local stress and strain asso- 
ciated with each cell. The cells elastically interact with 
each other; however, such interaction must be isotropic 
for the present analysis to apply. A configuration of dam- 
age is described as a damage field ip corresponding to the 
set of iV local variables (p = {ipx}x^n where N = Card(ri) 
is the total number of cells in the system. 

Our systems are initially uniform by hypothesis; i.e., 
they have a homogeneous damage field ip = at zero 
strain and stress and each cell starts with the same elas- 
tic moduli. Strain is progressively applied through the 
application of a uniform normal displacement I at the 
edges of the system. A cell breaks at constant appHed I 
when the energy required to break it (which is a random 
quenched threshold sampled from a probability distribu- 



tion function) just equals the reduction in stored elastic 
energy in the lattice due to the break. 

A key requirement of the class of models treated here 
is that the stress perturbation emanating from a broken 
cell must be weak enough that a first Born approximation 
holds. This means that the stress interaction between 
any two broken cells is allowed for while simultaneous in- 
teraction between three or more broken cells is not. This 
approximation is vaHd whenever the stress perturbation 
due to a broken site is much smaller than the mean stress 
in the system. 

One specific realization of such a "weak damage" model 
is an appropriately defined fiber bundle model. In the 
model, a set of TV elastic fibers are stretched between a 
free rigid plate and an elastic halfspace. The rigid plate is 
displaced by a controlled amount I that puts the fibers in 
tension. Once fibers begin to break, elastic interactions 
occur through the elastic solid. As demonstrated in the 
appendix, such interactions will be weak if either: (1) the 
elastic soHd is much stiffer than the fiber material; or, (2) 
the fibers are much longer than wide and are sufficiently 
widely placed. In the limit that the elastic halfspace be- 
comes rigid, this model reduces to the classic global-load 
sharing fiber bundle. 

Another reaHzation is a uniform elastic isotropic soHd 
divided into N cells. Uniform displacements I are ap- 
plied normally to the limiting faces of the lattice in such 
a way that the material is in a state of uniform dila- 
tion. The damage that arrives is assumed to change the 
isotropic moduH without creating anisotropy in the pro- 
cess. For example, the damage might be modeled as a 
spherical cavity that opens at the center of the cell thus 
reducing the elastic moduli of that cell. The weak inter- 
action is guaranteed if the change in the cell modulus is 
small. This can be considered a special case of the ori- 
ented damage models we have considered in earlier work 

As demonstrated in the appendix, the energy Ep re- 
versibly stored in such systems when the system is in a 
damage state (p at the applied loading level I is 

Ep[^J] = (Co + Ci+C2)f (1) 
Co ^ N (2) 

Ci = -c^(p^ (3) 

C2 = -e^Jxy(Px<Py (4) 

where c is a positive constant in the range < c < 1 that 
is independent of the damage/ deformation state, e is a 
small positive parameter in the range < e/c <C 1 that 
controls the strength of the stress perturbations, and Jxy 
are 0(1) coupHng constants that allow for the load re- 
distribution between cells at positions x and y when a 
cell is broken. Various spatial ranges for Jxy are con- 
sidered including: (1) global load sharing in which case 
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J'xy = 0; (2) local load sharing in which case Jxy = a > 
is a constant of order unity for pairs xy of nearest neigh- 
bors, otherwise; (3) elastic load sharing in which case 
Jxy ^ (^/k — y\)^ where i is the lattice step size; 
and (4) fiber-bundle elastic sharing with fibers interact- 
ing through an elastic plate in which case D — 2 and 
Jxy — y\ [see the appendix]. 

The local load sharing case (2) can in principle hap- 
pen in a fiber bundle stretched between plates enduring 
both elastic and plastic deformations capable of screening 
the stress perturbations caused by broken fibers to only 
nearest neighbors, which always can carry some load if 
damaging them only decreases their elastic constant (i.e. 
c < 1). We treat case (2) for the sake of generality and 
do not specify the detailed constitutive relations required 
for it to be realized in practice. 

The cost in surface energy to break any of the cells 
(which represents either the energy required to create 
new surface area within a cell or to break a fiber) is a 
random variable fixed ab initio, with no spatial correla- 
tions between the different cells. The breaking energy 
is thus a quenched uncorrelated disorder, described by a 
probability distribution function p{e) for which p{e)de is 
the probability that a cell's surface energy is in [e, e-|-de], 
and having a cumulative distribution P{e) = J^p{z)dz. 
For a given realization of each cell's surface energy Cx, 
there is thus a certain total surface energy 

Es[f] ^^ex{l ~ fx) (5) 

xen 

associated with each damage field f = {'Px}x£n- 

Given that B C fl is a certain subset of locations, 
the notation (p^ refers to the state where this subset is 
broken and its complementary intact; i.e., (px = 1 for 
every x ^ B, and (px = for every x G I = B'-^ . 

As the external deformation I is increased, damage evo- 
lution is ruled by Griffith's principle: Given that the sys- 
tem is in a certain configuration tp^ , at a given deforma- 
tion Z, it can undergo a transition towards a more broken 
state {p^^^^^ differing from the previous one by one ad- 
ditional broken cell at x, if the release in potential energy 
is equal to the surface energy cost of the new break; i.e., 
if 

/^Ep[p^,x,l]=ex (6) 

where 

Aii;p[^^,x,Z] = Ep[p^,l]^Ep[ip^'~'^-\l] 

= {c + eY,J^vVy)l''- (7) 

y 

If ex > AEp[ip^,x,l] for any surviving cell x £ I, there 
is no break and the deformation can be further increased 
while the system remains in the same state (p^ . If a break 
happens in cell x, which leads to the new state p^ where 



E = B U {x} while the external deformation / is kept 
constant, there is a possibility of avalanche at fixed I if 
there is some y € E^ such that 

AEp[p'',y,l]>ey. (8) 

If there is more than one possible location satisfying Eq. 
(jHJ, the one which breaks is determined by maximizing 
the energy release; i.e., its location corresponds to 

AEp[p^,y,l]-ey = M!ix,^EC {AEp[p^ , z,l] - e,) . (9) 

The avalanche test [Eq. 0] is then computed again until 
the system stabilizes in a given configuration. 

Eventually, although we have chosen to base the evo- 
lution of our damage model on minimization of energy, 
we note that the formal integration of path probabilities 
presented in the following sections could similarly be ob- 
tained as well for the case of a rule based on force thresh- 
olds, with at any given external deformation, a force car- 
ried per intact fiber equal to an average one, plus pertur- 
bations due to the already broken fibers. However, this 
approach will not be pursued here. 

III. PROBABILITY DISTRIBUTION OF 
DAMAGE STATES 

The probability of occurence of any damage configura- 
tion p is now determined when the system is at a given 
external deformation I that was reached monotonically {I 
as defined here does not include any elastic unloading). 
We will first summarize the results for the simplest case, 
global load sharing, which was performed in and 
which will serve as a basis for a perturbative treatment 
to include the effect of interactions. 

A. Global load sharing 

In this case, the interaction term of Eq. Q is C2 = 
for any configuration, and from Eq. 10, AEp[p,x,l] = 
cP regardless of the state p and new break location x 
considered. 

Each of the cells then share the same level of defor- 
mation I, and the probability for any one of them to be 
broken is simply (from the cumulative surface energy dis- 
tribution) 

Po P{e < cf) = / p{z)dz. (10) 
Jo 

The probability of being in a configuration p with n cells 
broken out of N is then 

p[pj] = Poiirii~Po{i)f-''. (11) 

This corresponds to the behavior of a so-called "two-state 
system" in which each of the TV independent sites has a 



probability Po{l) of being brok;en and 1 — Po(0 oi being 
intact. The consequences of this distribution function 
for the mechanical behavior and correlation length at the 
transition point of macroscopic rupture will be developed 
in section FV Al 



B. Local load sharing 

1. Damage-state probability distribution 

We now consider the case where each broken cell in- 
creases the local deformation by an amount ael on each 
of its nearest neighbors; i.e., Jxy = a for each nearest 
neighbor pair, or Jxy = for more distant cells. 

Note that for a given cell x, the potential energy release 
defined in Eq. JJjl, AEp[ip, x,l], is a growing function of 
both the deformation I and the subset of cracked cells ip 
considered - in the sense that if we consider two subsets 
B C D and a; G D'^, then AEp[ip^,x,l] > AEp[ip^ , x,l]. 
Physically, this inequality means that the local load over 
each cell increases with the external load imposed, and 
that any cell breaking anywhere else induces an addi- 
tional increase in local load. This inequality will play a 
key role in obtaining the damage-state probability distri- 
bution. 

To aid the pedagogic development, we first derive the 
probability of occurence of a configuration Lp^ with n 
isolated broken cehs forming a subset B, which do not 
share any common nearest neighbors. There are zn near- 
est neighbors to these n cells (in a subset F corresponding 
to the boundary of broken cells) where z is the coordi- 
nation number of the lattice considered. There are then 
N — {z + l)n cells completely isolated from the broken 
ones in the subset {B U F)'-^ . For a cell x having broken 
from an intermediate stage (p^ with ^ C B, at an inter- 
mediate external load ft. £ [0, ^], the change in the stored 
energy satisfies AEp[ip^, x, h] — ch? < cP . For every cell 
having survived, we have either AEp[p^ ,x, I] = {c-\-ae)P 
lix ^ F {x\s on the boundary of the broken cells set), or 
AEp[<p^,x, I] = cf ifxe{BU F)^ {x is completely dis- 
connected from the broken cells set). Applying Griffith's 
principle [Eq. Q] to every cell and intermediate defor- 
mations h g [0, /], and using the monotony of AEp[ip, x, I] 
in both I and p, we obtain that a necessary and sufHcient 
condition for ah cehs in B to be broken is that their sur- 
face energy thresholds were below d^, while those of their 
neighbors in F were above (c-|-ae)P, and the remaining 
ones m [B U F)^ above cP . Defining 



P^{1) = P{e < {c + ae)f) 



piz)dz, (12) 



this then imphes that 

P[(p, x] = P^'(l - Pi)^"(l - Po)^"'^+^^". 
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Figure 1: Typical configuration and associated connectivity 
measures. 



is the probability of occurence of such a configuration. 

We now pass to the more general case. In the argu- 
ment, we obtain upper and lower bounds for the proba- 
bility of some arbitrary damage state, and then demon- 
strate that in the limit e/c <C 1, the two bounds converge 
to the unique probability distribution of interest. 

For any configuration ip, Uk is defined as the number 
of intact cells with k broken nearest neighbors and nk 
as the number of broken cells with k broken neighbors. 
If n out of the N cells are broken, we have J2k ~ ^ 
and Uk ^ N — n. The way the above quantities are 
associated to any particular configuration is illustrated 
in Fig. El 

For a cell x which broke from an intermediate state C 
p with k already broken neighbors in ■(/'; at an external 
load h, we have 

AEp[il:, X, h] = {c + kae)h^. (14) 

Using again the monotony of AEp, a necessary and suffi- 
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cient condition for any intact cell y € B'^ to have survived 
is that its threshold exceeded Ai?p[iy9, /]. The proba- 
bility for each of these independent statistical events to 
occur is expressed as 1 — Pk{l), where 



Pk{V) = P{e < {c + kae)l^). 



(15) 



For any cell which broke x e B, we note that they have 
broken with certainty at the ultimate deformation / if the 
external load was sufficient to trigger their break without 
the help of overload due to breaks of the other ones; i.e., 
they have broken with certainty if their energy threshold 
was below AEp{ip'^, x, I) — cP where denotes the empty 
set (no broken cells). Furthermore, if every threshold in 
B was below cP, except a particular one x £ B which has 
kx broken neighbors in the considered configuration ip^ 
and has a breaking energy between cP and (c + kxae)P, 
the kx neighbors of this considered cell x will have bro- 
ken with certainty at the ultimate load I, so that x will 
also break with certainty under the effect of the overload 
due to its broken neighbors. The probability that this 
individual threshold broke under the sole effect of the 
overload due to its neighbors can be expressed 



PkAi)-Po{i) 



(16) 



Thus, a lower bound for the probability of occurence of 
the configuration (p^ can be expressed as 




Uil-Pk)'"' (17) 



keN 



where the index k runs formally to N; however, — 
when k > z (the coordination number of the lattice). 

For any cell x & B which broke, its associated energy 
threshold was necessarily lower than AEp[(p^'^^^\x,l] 
where B \ {x} denotes the set B with cell {x} excluded 
from it. Thus, an upper bound for the probability of 
occurence of the configuration under study, is 
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11(1 - P,)-" 



(18) 



For a continuous p.d.f. over thresholds (that guarantees 
no jumps in the cumulative distribution P and is the 
only restriction placed on the p.d.f.), APfc is a quantity 
of order e, and Pq is of order 1. In the limit of a Born 
model (weak stress perturbations for which e/c ^ 1), the 
upper and lower bounds for the probability of occurence 
are identical to order e. We have thus established in this 
framework that 



P[ip,l]^ l[Pm''-l[{l-Pkr^ 



(19) 



m=0 



k=0 



2. Identification of a surface tension and cohesion energy 

The above can be re-expressed for small interactions 
e ^ 1 by a Taylor expansion of the cumulative q.d. 
(quenched disorder) distribution as 



APfc(/) = 7(0fc where 7 = p{cP) aef . 
Then, Eq. ifTfljl becomes 
lnP[^,/] - 
In 



(20) 



X{{p^+imr-X{{i- 

.m=0 fe=0 

nlnPo + (iV-n) ln(l - Pq) 



Po - 7^)"'' 



(21) 
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nlnPo + (iV-n) ln(l 
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" fc=i 
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2—ni 
Po 



1 -Pr 
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(22) 



where in the state ip considered, n is the number of broken 
sites, and ni and ns refer respectively to the number 
of internal bonds between pairs of broken cells and the 
number of boundary bonds between broken and intact 
cells. 

These probabilities are thus of the form 

P[ip,l] = p[(^®]e"'''(^°/(^--f'o))+2"^'>'/-Po-"s7/(i-fo) (23) 

with the probability of the intact state 

P[p^] = (1 - Po)^. (24) 

The probabilities are thus classic Boltzmann distribu- 
tions P[tp'^]e^-^ where 



-H[^,l]^ In 



Po 



(1-Po 



2—ni 

Po 1 



Pf 



■ns- (25) 



For no interactions, 7 = and the above reduces to the 
global-load model of Eq. lfTT|l . Since n and nj can be 
formally interpreted as volume integrals over the interior 
of clusters of broken cells, and ns as a surface integral 
along the boundary of these clusters, we can make the 
following analogies to the quantities of classical statistical 
physics: 



In [Po/(l-Po)] 

^0 

7^ 



-^J^|T chemical potential (26) 



1-Pn 



T 
ls_ 
T 



bulk energy (27) 
surface tension. (28) 



The first term in Eq. lf25|l accounts for the average energy 
required to break a cell, the second term for an increase 
in the probability of finding some clusters of connected 
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cracks due to positive interactions between them, and 
the third term for a decrease in the probabihty of finding 
clusters with a long interface between cracked and non- 
cracked region due to the fact that intact cells along the 
boundary are more likely to have broken from overload- 
ing from cracked neighbors, thus leading to even more 
fractured states. 

In the beginning of the process, the two first terms 
dominate, and to leading order 



-H[^,l] = {lnPo)n + 2^ni. 

-TO 



(29) 



For example, defining S — limj^^o[lnP(2/)/ln2/], so that 
P{y) ^ ay^ , in the small deformation Hmit, one has: 

a5ec~^\ and, 
ln(ac*P'^)ri + 



Po ^ ac^r"; 7 ~ aaSec^'-H''''; 7/P0 ' 
in decreasing order of magnitude, 



ac^P^n + 2a5ec- 



aaSec 



ns. 



C. Arbitrary- range interactions 

As long as screening effects are absent or neglected (as 
they are in the present model of weak interactions) , the 
above arguments based on Griffith's principle and the 
monotony of AEp in f and I, extend directly to the case 
of arbitrary ranges of interactions. 

At external deformation I, a cell x has: (1) broken with 
certainty if the external load alone could break it; i.e., if 
its associated surface energy is lower than cP; and (2) 
is intact with certainty if the external load plus the load 
perturbation due to the broken cells in the final configu- 
ration could not break it at final deformation; i.e., if its 
surface energy is higher than {c + eJ2y Jxy'Py)^- Denot- 
ing 



Po 



P{e < cP) (30) 
P(e < (c + e ^ JxyVy)l^) - P{e < cf) 



(31) 



a necessary condition to end up at a certain configu- 
ration Lp^ at deformation I is that all surviving cells 
X e B'-^ in that configuration have their threshold above 
{c + eY^y JxyiPy)P, and all cells which broke x £ B have 
their threshold below (c + £j2y Jxyfy )P- Thus, 

nv^'A <X{{Po + APx) n (1 - ^0 - ^^-) (32) 



xeB 



provides an upperbound for the probabilities in the case 
of arbitrary ranges of interaction. 

Conversely, a sufficient condition to end up in this con- 
figuration is that all surviving cells have their threshold 
above {c + sj^y Jxy'Py)P , and that the broken cells ei- 
ther have: (1) all their thresholds below cP; or (2) all 



but one located at x have such thresholds, the last one 
breaking only due to the overload from the others; i.e., 
the last one has its threshold bounded by cP < e < 
(c -I- Jxy'Py)P- This gives a lower bound for the 
probability of the configuration (p^ , 

Pb^,Z] > ^(Po+AP,) n Po Y[il-Po-AP,). 

xeB yeB^{x} zeB^ 

(33) 

As earlier, both lower and upper bound coincide to 
order e, so that upon Taylor expanding hiP[tp,l] to 
this order, we again obtain the Boltzmannian P[p>, I] — 
P e^-^I"^''! with an intact probability given again by 
P[/] = (1 - Po)^ and 

-H[ip,l] = ln[Po/(l - Po)]n 
picP)eP 



Po 



p{cP)eP 
1-Po 



^ ^ Jxy^xPy 
xy 

■'^Jxy'Px{l~ iPy)- (34) 



In the beginning of the process, the two first terms once 
again dominate and 



-H[p^,l] = ln[Po]7i 

p{cP)eP J2xy Jxy<PxVy 
-\ 



(35) 



This expresses the equivalence between this most gen- 
eral weakly interacting damage model and an Ising model 
with generalized interaction rules. 



IV. EQUIVALENCE WITH A MAXIMUM 
ENTROPY POSTULATE 



We now obtain these same probability distributions 
using the standard entropy maximization argument. 

It is convenient to introduce the index j to denote 
each possible damage configuration ip. We postulate that 
the probability distribution function over configurations 
j maximizes Shannon's disorder 



S ^ -^pj lnpj 



(36) 



subject to the constraints 

Y^PjEj = U ■ ^PjTij = n; Y^pj = 1; Vj Ij = I (37) 

j j 

where U is the total average energy that has been put 
into the system and n is again the average number of 
broken cells. The vaHdity of this maximization postulate 
will be directly demonstrated in what follows. However, 
independent of the formal demonstration, one can antic- 
ipate that Shannon entropy should be maximized since 
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the initial quenched disorder in the breaking energies al- 
lows each possible damage configuration to be accessible. 
The constraints allow for the content of the Griffith prin- 
ciple and are what make certain emergent damage states 
more probable than others. 

Throughout the remainder of the paper, pj denotes 
the probability of finding a configuration j over all pos- 
sible realizations of the q.d., when the system has been 
brought to average deformation I starting from an ini- 
tially intact state. The term rij now refers to the total 
number of cracks in the configuration j, while n is the 
statistical average of rij. A few classical results can di- 
rectly be derived from these assumptions: 



where translational invariance of Jx 



</r=v-x has been 



-fi{Ej~imj) 



z 



where P = dS/dU\i . 
these one further has 



z 



(38) 



and = — dS/dn\i ^. From 



G = -Hz)/p 

G = U-S/P-^ln 

dG = Tdl - Sd{l/I3) - ndn 

dU = Tdl + dS/P + ndn 



where 



Tj and Tj 



dE^ 
dl 



(39) 



(40) 



assumed. Equation Il4ip then requires 



P = 



M = 



Foil 
Foil 



Po) 
Po) 



picP) 



In 



(l-^o) 



(43) 



Thus, the p.d.f. over configurations is indeed maximizes 
entropy [Eq. Ij36|l ] under the constraints of Eq. ()37ll . with 
no unknowns. The inverse temperature and chemical po- 
tential depend on the deformation level through Eqs. 
lll3ll - l(iijl . They are well defined analytical functions 
of I and the model parameters considered, the q.d. dis- 
tribution via Pq, and the interaction coupling Jr. So 
the usual machinery of equilibrium statistical mechanics 
[Eqs. Ij38ll - lj39|l ]. is valid and can be used for any of our 
damage models. 

The autocorrelation function {ipxfy) can therefore be 
obtained in the standard way by: (1) defining a new 
Hamiltonian = Ej + QxVx that incorporates a 
coupHng with a formal external field Qx\ (2) evaluating 
the associated generaHzed partition function Z'; and (3) 
performing second-order derivatives with respect to the 
external field, 



1 



d^z' 



p-^Z' dQxdQy 



(45) 



The thermodynamic parameters P and /i are obtained 
here by comparing Eq. Ij38|l to the earlier probabilities 
obtained by direct integration over the microstate space. 

The p.d.f. over configurations is a maximum of Shan- 
non's entropy under the above constraints, if and only if 
there are two constants (/3, p) such that 



Vj, Hj — P{Ej — jirij) = constant 



(41) 



with Hj given by Eq. l|34)l . From Eqs. (QJ-Q and Eq. 
lIMji . we have 



1 

— ^)%- 

p{cf)eP 



E, 



^0 

xy 

p{cP)eP ^ 

^ _ p 2^ .JxyVx{l - fy) 

{N ~ CTlj - Jxy(pxiPy)l^ 
xy 



'^Jxy(Px{^~fy) = -^Jxy^fix^Py + rij'^Jr (42) 



xy 



xy 



Second order derivatives of the free energy G with respect 
to P and fi can also be directly related to variances of 
the number of broken cracks nj and the potential elastic 
energy Ej, but these standard derivations will be left to 
the attention of the reader. 

Note that the formal temperature (the energy 
yardstick used to distinguish probable from improbable 
states) behaves regularly throughout the damage process 
(for continuous q.d. distributions). In the case of a uni- 
form q.d. on [0,1], it will take the particularly simple form 
l/P = Po(l — -Po), starting from zero and going back to it, 
with a maximum when Po = 1/2. The chemical potential 
/i behaves regularly as well. 



V. APPLICATIONS 
A. Global load model 

A consistency check of the results in Section ITVl is now 
performed for the case of the simple global-load sharing 
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model. From Eqs. with Jr = 0, we have 



P{Ej — iirLj) = —In 



{N-cnj)l^ 



(46) 



1-^0 



hi 



Po{l - Po) 

Po 



l-Po 



[{N - cnj)P + cujl^] 
n,+f{l). (47) 



Independently, we also have the direct result from Eq. 



- Po"^(l - Po)"^""^' = Poe^"^'"[^°/(^"^°)l. (48) 



Thus, the Boltzmann distribution Eq. Ij^i8|l with temper- 
ature and chemical potential given by Eqs. (O and lO 
is indeed identical to the known solution Eq. lli8)l . which 
confirms the validity of the expressions for (3 and fi. 

It is also instructive to look at all terms in the first law 
Eq. l(39jl to see exactly what they represent. The values 
of the average quantities in this simplest model can be 
obtained using the lemma 



5]Po"^(l-Po)^-"^n,=7VPo 



which is demonstrated by applying the operator xd/dx 
to the binomial theorem 



{x + y) 



N 

N \ " 



^ [N -n)\n\ 



and then taking x = Pq and y = 1 — Pq. Using Eq. dV A.\ . 
one obtains 



n 
U 



PqN 

^Po"^(i-Po; 



N-r. 



N{1 - cPo)f 

-iV[PolnPo + (1 - Po) ln(l - Pq)] 



^ = {N-cn,)2l 
dl 

N{1 - cPo)2l. 



Taking the derivatives of these quantities yields 



dn = NdP(t 
dS 



(49) 



In 



Po 



fidn 



lidn - 
dS 



fidn 
dS_ 

T 

dS 

Tdl 

\- Tdl 



l-Po 
Pojl-Po) 
p{cP) 

Po{l~Po) 



dPo 



Po 



l-Po 
Po 



l-Po 



p{cP) 
= -cfNdPo 
= N{l-cPo)2ldl 
= N{1 - cPo)2ldl - cfNdP(_ 
= d[N{l - cPo)?2] = dU 



~ cPj NdPo 
NdPo 



which is a consistency check for the validity of the first 
law [Eq. The average mechanical behavior of this 

model, as well as the evolution of entropy and formal 
temperature are illustrated in Fig. |2l for fiat q.d. distri- 
butions between and ^max and for various values of c 
(the parameter that controls the relative change in stiff- 
ness due to a cell breaking). 

Most importantly for our present purposes, since the 
probabity Pq of having a site broken in this model is 
independent of the configuration and site location, the 
global load-sharing damage model is exactly equivalent to 
the percolation model with occupation probability Po(0- 
There is a critical-point phase transition in this model, 
when Po(^c) — 1/2, for which S goes through a maximum 
5" = iVln2. The correlation length diverges then as 



^^\Po{i)-i/2r 



(50) 



with ly 
general 



and thus 



4/3 in dimension D — 2 We have in 



Po{l) - Poilc) - P{cll){l - Ic) 



i-\l-lc\-'/'. 



(51) 



(52) 



In pathological cases, special q.d. distributions satisfy 
p{dl) = so that Po(0 - Po(Zc) ~ - with a I. 
This results in 



^-\l-lc\ 



-4a/3 



(53) 



Note that in Ref.jllJ, we have treated this model with 
c = 1, which did not change anything to the nature of 
the transition, but only changed the minimum stiffness 
associated with the most damaged configuration where 
Tmin — 2N{1 — c)l, and, consequently, the position or 
existence of a peak stress in the average mechanical re- 
sponse t{1). This is seen in Fig.^tb): the existence, and 
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critical-point nature of the percolation transition are in- 
sensitive to c. 

The present approach does not allow a direct explo- 
ration of the avalanche distributions as the critical point 
is approached. This is because the probability distribu- 
tion over configurations was obtained at each elongation 
I by averaging over all realizations of the quenched disor- 
der. To obtain directly a result on avalanches, on the con- 
trary, correlations between successive elongations should 
be considered for each realization of the q.d., and the 
average over q.d. should only be considered afterwards. 
Our damage model nonetheless behaves as a fiber bundle 
with GLS for which results about the avalanche distribu- 
tion have been determined (hIIi^IT^ . If c is sufficiently 
large that the average strain/stress curve has a peak at 
Tc, there is a burst in the number of elements that must 
break before the system is able to sustain the same global 
force again. The distribution of the burst size A should 
thus scale as 

n(A) ~ A-5/2e-A/Ao^ 

with Ao ^ |r — rd^^/^. It might be tempting to re- 
late this power-law behavior of the burst distribution to 
the behavior of cluster numbers in mean-field percolation 
(on a Bethe lattice or on hypercubes in dimensionality 
-D > 6). The average number of clusters of size s per 
lattice site also scales as rip (A) ~ lsr^/'^e~^^ ^° in such 
models where Aq ^ \p ^ Pc\^'^ i^k- -^^t the fact that 
both the cluster number distribution at critical perco- 
lation threshold and the burst size distribution at peak 
stress have the same power law exponent of -5/2 might 
simply be coincidental. These two points (critical per- 
colation threshold and possible peak stress) in general 
differ, and the exponents characterizing the evolution of 
the upper cutoff Aq of these distributions close to perco- 
lation threshold or peak stress seem also unrelated. 



Local load model 



Figure 2: Thermodynamic and mechanical response as a 
function of imposed deformation for the global-load sharing 
model: (a) the difference between the average Young's modu- 
lus T /I and the minimum Young's modulus Tmin/l (that holds 
in the entirely damaged configuration); (b) average stress for a 
few particular models (the black dot is the percolation critical- 
point transition); (c) Shannon entropy; and (d) temperature. 



In this case, Jxy — a when x, y are nearest neighbors, 
and zero otherwise. Defining the order parameter a G 
{ — 1, 1} as cr = 2(/j — 1, we note that 

'^Jxy'-Px'-Py ^ OiN + 2a'^ax + ^ <yx<yy (54) 



xy 



<xy> 



the position of a possible peak-stress position relative to 
the percolation transition (black dot l/lmax ~ 1/v^)j de- 
pends on the particular model considered. However, the 
divergence of the correlation length and the associated 



where < xy > denotes a sum on nearest neighbors only. 
Consequently, the probability over configurations can be 
cast under the form 



■,P{ET.:c'^^+JT.<^y>'^^'^y) jZ^ (55) 
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which is exactly a classical Ising model with coupling 
constant and external field given by 



f3J = 



(3E 



4Po{l-Po)' 



1 



In 



1-Pn 



2ap{cP)ef 



(56) 
(57) 



The critical point of this model is at ;64| 

(/3J„/3i?) = (A„0) 

with Ac = ln(l + V2)/2. The external field PE starts at 
infinitely negative values, and ends up at infinitely pos- 
itive ones. It evolves continuously and thus necessarily 
crosses i? = at the Ic satisfying , from Eq. 115 711 . 

Poilc) {ln[l - Po{lc)] - ln[Fo(^c)]} = 4.ap{cll)ell (58) 

The mean-field percolation result of Po(^c) — 1/2 is thus 
recovered when the coupling vanishes (J — 0), which is 
a consistency check. 

More generally, for non-zero nearest coupling constants 
a, the system will undergo a first-order transition if at Ic 
satisfying Eq. Ij58ll , the formal inverse temperature satis- 
fies 



l3{lc)J 



ap{dl)ell 



4Poilc)il - Poilc)) 
ln[l-Po(g]-ln[Po(/,)] 

16(l-Po(/c)) 



> Ac 



(59) 



Depending on the value oi f}{lc)J, the system can display 
four types of behavior, that are schematically depicted in 
Fig. El 



t»3 



(Hi) 



(ii) 



\ (iv) 



Ising 
critical 
point ^ 



First order transition line 



Pi 



Figure 3: Possible Paths in the space of coupling constants of 
Ising models, under increasing imposed elongation associated 
to different Local Load Sharing damage models. 



(i) For many q.d. distributions, the first value of Ic sat- 
isfying Eq. (I58|l occurs for very small values of Poih), 
that correspond to small values of 2ap(cZ^)eZ^, since 
xln[(l ~ x)/x] — > O"*" when a; ^ 0. In this case, 

in[(i - Po(;,))/Po(/c)]/(i - Po(U) - -HMQ] > 1, 

and the condition of Eq. Ij59|l is fulfilled. The system 
thus goes through a first-order phase transition at this Ic 
and there is a discontinuous jump in the average number 
of broken cells and the average stress (which are related 
to first derivatives of the free energy with respect to I or 
/i, and are similar to the average number of spins up in 
an Ising model) [3| . The correlation length increases up 
to the transition but remains finite. All of this behavior 
for such local load models has been documented in the 
literature [HElEi- 

However, three other behaviors are also possible if 
P{lc)J < Ac at the point Ic at which the external field 
vanishes [Eq. 15811 ]. Whether this occurs is controlled by 
the type of q.d. distribution and the value of the coupHng 
constant a. When (3{lc)J < Ac, no first order transition 
is encountered and one can further classify the transition 
into three subcategories: 

(ii) If f3{lc) has a finite value of order unity, signif- 
icantly below Ac, the system simply goes continuously 
through Ic, without discontinuity in sustained load or av- 
erage number of cracks. The correlation length remains 
finite. 

(iii) If l3{lc)J ^ 1, (which should happen for vanishing 
a), the distribution over configurations is dominated by 
the external field, and the system essentially behaves as 
a percolation model going through the percolation tran- 
sition. Although there should be corrections due to the 
nonzero character of (3(lc)J , these might be smaller than 
the finite-size corrections in numerical realizations, and 
the correlation length would then be found to diverge up 
to the system size as ^ ^ |Z — lc\~^^^- 

(iv) Finally, \i (3{lc)J A' , the system comes close to 
the critical point of the Ising model at Ic, and the corre- 
lation length should diverge corresponding as in an Ising 
system whose temperature comes close to 1/ Ac — Tc from 
above, while the external field reverts sign. The slope 
of the average mechanical curve t{1) should also locally 
diverge around Ic- The exponents associated with the di- 
vergence of ^ as function of {I — Ic) depend on the way the 
critical point is approached as a function of I. For values 
of I such that (3E{1) < 1, we write l/P{l)J ^Tc + f{l), 
and C - |r(/) - Tc]-" = 1/(01"'' where v = \ ior the 2D 
Ising model. The correlation length therefore diverges as 
\l — lc\~'' , unless the temperature has a quadratic mini- 
mum in I close to Ic in which case it diverges as \l — lc\~'^'^. 



C. Power- law decay 

We last consider the general case of stress perturba- 
tions decaying as power-laws of the distance to broken 
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cells, J{x) a{x/d)~"' /^(j), where a is of order unity, 
d is the lattice constant, and $(7) = J2y^x(\y ~ ^\/^)~"' 
is a normalizing factor (which depends on the lattice size 
if 7 > D, where D is the system 's dimension). To be con- 
sistent, the models considered here are biperiodic of lin- 
ear size L, and the interactions are put to for distances 
above L. This type of model, considered for example by 
Hidalgo et al |lj,lla|, allows us to span ranges between 
purely global sharing (when 7 = 0), to the local sharing 
limit (7 — > 00). Equations H42I44|I show that this model 
leads to probability distributions over configurations of 
the form 



P[a,l] 



/Z, (60) 



which is a generalized long-range Ising model with cou- 
pling constants and external field given by: 



/3J(r) 



4Poil-Po) Hi) ' 



/3E = -ln( — 



2Fn 



(61) 
(62) 



Note that although the case 7 ^ 00 is isomorphic to the 
Local Load Sharing Model introduced in the previous 
section, the Global Load Sharing Model of Section FV Al 
corresponds to a = 0, but not to the exponent 7 = 0. 
The presence of the quadratic coupling makes this later 
case isomorphic to a Curie-Weiss model, which is the 
Mean-Field theory of Ising models. 

Although such long-range Ising models are still an open 
area of research, it has been proposed that 
their behavior can be classified into two categories de- 
pending on whether 7 > Z) or 7 < Z?. In the short- 
range case, J > D, this model admits a traditional ther- 
modynamic limit when L 00. Although $(7) di- 
verges when L 00, the thermodynamic limit is well 
defined once $(7) is introduced into the coupling con- 
stants J(r) = (r/d) (7). Thus, the free or internal 
energy per lattice site, entropy and magnetization all ad- 
mit a finite limit when L ^ 00, and are functions of /3 
and the external field E. 

In our case, we have 



/3 
(3E 



ap{cP)eP 
4Po(l - ^0) 

iln(^' 
2 ^ 1 - Po ' 



ap{cP)ef 
2Po 



With J defined as above, the Ising model shows that in 
the limit 7 ^ 00, the critical point is at 



lim Pc 

7 — >-oo 



ln(l 



2)/2. 



It has also been shown that as 7 ^ D+, f3c ^ I- 
All short-range models 7 > D display critical points 
at finite temperature, and since all of them have a fi- 
nite range of interaction, one can conjecture that all of 



them belong to the same universality class as the Ising 
model, with a critical point such that Pc is of order 
unity [iJl. It has also been conjectured that all non- 
local models 7 < D belong to the Curie- Weiss uni- 
versality class of the Mean Field Ising model. They 
also present critical temperatures such that Pc is of or- 
der one The Curie- Weiss model corresponds to 
the mean-field coupling Jxy = ^/N, independently of x 
and y, where N is total number of cells in the system. 
For the Curie- Weiss model, the p.d.f. over state config- 
urations is P[a] OC exp[/3(£'X]2; O-a; + J2xy'^xO-y/N)] cx 

exp[/3(i?E.^. + (E.<^.)V2A^)]- 

All this suggests the following classifications of our 
damage models: 

(1) If the coupling constant is small (a <C 1), the dam- 
age model is close to the percolation model so that in the 
approach to the transition at elongation Ic, the correla- 
tion length behaves as ^ ~ |/ — lc\~'^ , with v = 4/3; 

(2) For non-negligeable coupling constants a, in the 
short-range case j > D, we recover the same three possi- 
ble behaviors as described above for the local load sharing 
rule; 

(3) In the long-range case where 7 < Z?, we again re- 
cover the same type of four scenarios depending on the ra- 
tio P{lc)/Ac (as discussed in the local-load sharing case). 
The model behaves as percolation when a ^ 1. Other- 
wise, the behavior is determined by the ratio P{lc)/Ac. 
Note that Ac is of order unity, but depends on the partic- 
ular exponent 7 and on the system size. If P{lc)/Ac < 1, 
the system behaves continuously and no transition is ob- 
served. If P{lc)/Ac > 1, there is a discrete jump in both 
the average number of broken cells and the sustained 
load. The correlation length remains finite in both of 
these two cases. Only the limiting case of P{lc)/Ac ~ 1 
corresponds to a second-order phase transition. In any 
of these cases, for large enough systems, the Curie- Weiss 
description holds according to Refs. [63,^fi^E3- Accord- 
ingly, if there is any divergence of correlation length due 
to the system coming close to P{lc)/Ac ~ 1 at i5 = 0, the 
associated exponents should be those of the Curie- Weiss 
mean-field critical point {v = 1/2), and not the Ising one. 

These results can be compared to numerical-simulation 
results in the literature for related models. In fiber bun- 
dle models with power-law interactions , a transition 
has been found as a function of the interaction exponent 
7 that is consistent with the above analysis, predicting 
mean-field behavior for the long-range case 7 < Z?, and 
Ising-like behavior in the short-range case 7 > Z?. Typ- 
ical configurations prior to breakdown for this type of 
system are displayed in Fig. 5 of Ref. and look very 
similar to percolation configurations close to the tran- 
sition in the case 7 = 0, displaying shorter and shorter 
cluster sizes (characteristic of the autocorrelation length) 
and compact configurations as 7 increases above D. This 
is coherent with a mean-field behavior close to percola- 
tion transition in the first opposed to a first order 
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transition (analogous to Ising model crossing E = be- 
low Tc) when j < D. This analogy is even more apparent 
in Fig. 7, bottom, of Ref. where an extension of this 
model was considered, with time-delayed fiber breaking 
process in addition to this power-law decaying interac- 
tions [lg|- 

Burned fuse models, in which the interactions between 
burned fuses decay as , exhibit diverging autocorre- 
lation lengths at breakdown, with ^ ~ |^ — ^d"" where 
I is the imposed voltage, and ly is equal to the percola- 
tion exponent 0, v = 4/3 in 2D Q or 0.88 in 3D [s3|. 
The morphology of the connected "fracture" at break- 
down is oriented, and different from percolating clusters 
in the percolation model. This seems to be related to the 
anisotropic character of interactions in the burned-fuse 
model; i.e., the current perturbation from a burned fuse 
goes as a dipolar field decaying as l/r^, but also hav- 
ing an orientational aspect not included in the models 
under study in this paper. Based on this, Hansen et al. 
[i^ relate the roughness of the spanning fracture to the 
autocorrelation length exponent, based on arguments of 
percolation in gradient, and this properly predicts the 
Hurst exponent of the final damage fronts, both in 2 and 
3D. This anisotropic aspect is absent from the models dis- 
cussed here, but the fact that the autocorrelation length 
exponent is similar to the ones of percolation, is coherent 
with the fact that long-range systems are in the perco- 
lation universality class. Burned fuse systems are at the 
verge between short and long range interactions, in the 
sense that they correspond to -f = D. 

Last, fiber bundles connected to elastic plates, where 
7 = 1 and D — 2, have been numerically studied in Ref. 
[231 and an autocorrelation length exponent of ~ 1.54 
numerically determined. The present theory does not ex- 
plain this autocorrelation exponent. The discrepancy be- 
tween this result and the percolation or Curie- Weiss crit- 
ical point result, might presumably result from finite-size 
effects making this model (7=1) still significantly dif- 
ferent from the Curie Weiss one (7 = 0), or from the fact 
that the stress perturbation was in this numerical model 
too large for the small perturbation expansion performed 
here to apply. 

VI. CONCLUSIONS 

We have treated a class of damage models having weak 
isotropic interactions between cells that become damaged 
in the lattice. A quenched disorder is present in the rup- 
ture energies for each lattice cell, and the evolution of 
damage is ruled by the Griffith principle. Averaging over 
all possible realizations of the underlying quenched disor- 
der, the probability distribution of each possible damage 
configuration was obtained as a function of the deforma- 
tion externally applied to the system. 

The exact calculation is analytically tractable in the 



case of a global load sharing model, and it has been 
shown to be isomorphic to a percolation model. This 
corresponds to the behavior of a system totally domi- 
nated by the underlying disorder, where the next cell 
to break is always the weakest one. Spatial interactions 
added to the system modify this picture, since the over- 
load created by broken cells induce some degree of spatial 
ordering that competes with the weakest-cell mechanism. 
By limiting ourselves to small overloads compared to the 
average load of the system, it was possible to obtain the 
probability of damage configurations as integrated over 
all realizations of the quenched disorder. 

In this weak interaction limit, the resulting proba- 
bility distributions were shown to be Boltzmannians in 
the number of broken cells and in the stored elastic en- 
ergy. This type of distribution maximizes Shannon's 
entropy under constraints related to the energetic bal- 
ance of the fracture process, and we have demonstrated 
the formal relationship between our quenched-disorder 
damage models and the standard distributions arising in 
equilibrium statistical mechanics. This then allows the 
standard toolbox of statistical mechanics to be applied 
to our damage models. 

Our systems map onto three types of possible behav- 
iors: (1) percolation models in the case of interactions 
so weak, they may be neglected; (2) Ising models for 
non-negligeable short-range interactions; and (3) Curie- 
Weiss mean-field theory for non-negligeable long-range 
interactions. The temperature and external field in the 
partition function of our models are analytical functions 
that depend on the particular sharing rule, on the type of 
quenched disorder considered, and on the average elonga- 
tion (or deformation) externally loaded onto the system. 
The path followed in the Ising control parameter space 
when the load is increased from depends on the q.d. 
distribution and the load-sharing rule. When the formal 
external field reverts sign, a phase transition is possible. 
This can correspond to a first-order phase transition, a 
percolation-like transition, or an Ising critical-point tran- 
sition, depending on the value of the formal temperature 
during the transition. 

The systems studied here are limited to isotropic load 
perturbations. We have earlier studied oriented crack 
models in [i3,I^E3|' which corresponds to anisotropic 
load perturbations that depend on the orinetation of 
the crack opened in the lattice. Those earlier studies 
were based on an entropy-maximum assumption. The 
hypotheses of the present work extend directly to ori- 
ented systems, and so the present paper justifies the 
entropy-maximum assumption postulated in our earlier 
work. The precise value of the temperature, and the 
physical interpretation of the functional forms given in 
[Hnl I51I IH^. should be modified according to the results 
of this paper. Such modifications will, however, result in 
identical functional forms relating the configuration space 
and the p.d.f. over configurations, and thus the present 
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work confirms the existence of a phase transition in such 
an oriented crack model, with an associated divergence 
exponent of the autocorrelation length, v — 2. 



Appendix A: RECOVERABLE ENERGY AS A 
FUNCTION OF THE DAMAGE STATE 



The argument here will be specific to a fiber bundle 
model. However, as noted in the text, other weak damage 
models will also be controlled by the same type of stored- 
energy function obtained here. 

Define a fiber bundle as (initially) TV fibers stretched 
between a free rigid plate and an elastic half-space. The 
rigid plate has a controlled displacement / applied to it 
that stretches the fibers and the elastic half-space. As a 
fiber breaks at fixed I, the force it held will be transmitted 
to the other fibers through the elastic half-space. 

A fiber at point x is stretched a distance i^. Where 
that fiber is attached to the elastic halfspace, the surface 
of the halfspace displaces an amount Ux- Thus, at those 
places X where fibers exist, we have 



I 



(Al) 



The fiber at point x exerts a force on the elastic half- 
space that is 



Lp 



(A2) 



where A/ is the x-sectional area of the fiber (assumed to 
be independent of the extension), Li? is the inital length 
of each fiber, and Yp is the Young's modulus of each 
fiber. The local order parameter (^x is if the fiber is 
intact and 1 if broken. 

The Green function for point forces acting on the sur- 
face of an elastic halfspace ^69] yields 



W3(x) 
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where Y^ is the Young's modulus and the Poisson's 
ratio of the elastic solid. In general, the displacement at 
a point x = {xi,X2,x^) within the elastic solid (where 
— Q defines the surface) due to a point force acting 
at a point y = (yi,?;2,0) on the surface [i.e., f(x) = 
/y'^(x- y)3] is given by 
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where 



r = [{xi - yi)^ + {x2 - y2f + xl 



211/2 



(A8) 



Putting a;3 = in the expression for U3(x;y) and then 
summing over all y yields the expression for the displace- 
ment Mx of the surface. 

We now define the dimensionless number 



^ _ (1 - al) Yp Ap 

TT Ys LpLp' 

where a length Lp has been defined as 



1 



E 
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i.e., this sum is independent of which point x is consid- 
ered. Assuming either that the elastic halfspace is stiffer 
than the fibers, or that each fiber has a length much great 
than its width, or that fibers are spaced far enough apart 
that Lp is large, allows e to be considered a small num- 
ber. Since the fiber bundle is assumed to be made of a 
finite number N of fibers, there is no divergence to L_p. 

Using these definitions along with Ux = I — £x and 
iterating Eq. IIA4|I once to get the leading order in e con- 
tribution gives 



|y -X 



■^y +0(e2). (All) 



The elastic strain energy reversibly stored in each surviv- 
ing fiber is then 



1 ApYp 
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where terms of O(e^) have been dropped. Thus, upon 
summing over all the fibers we obtain the total energy 
Ep stored in the fibers as a function of the damage state 



(1 - 2e)N - (1 - 4e)^(^x - ej] J^yV^Vy 

X x.y^^x 



where the coupling constant is defined 
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This form of the fiber energy is consistent with what was 
defined in the text. 

We now demonstrate that the energy recoverably 
stored in the elastic solid makes no important modifi- 
cation to Ep. The strain energy in the solid is given by 



2(l + a,)l3>o 



i — la. 



(A16) 

where summation over the indices is assumed and where 
the strain tensor is defined 



eifc(x;y) = 



1 / aui(x;y) 9ufc(x;y) 



2 V dxk 



dxi 
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The displacements are given by Eqs. ifl^ll-lfTzIl . 

Prom these equations, the strain at points x inside the 
elastic soHd takes the leading-order in e form 



dk — elLp 



(A18) 



where the constant tensor Cjfc has units of inverse-length 
squared and the average strain tensor throughout the 
elastic solid is elLpdk. The perturbation term due to 
broken fibers volume integrates to zero. The tensor /it- 
has no dependence on the norm |x — y|; however, this 
fact is immaterial since fik plays no important role. 

Upon forming the required products for the integrand 
m Eq. PT6|l . and using the fact that terms linear in the 
broken-fiber perturbations integrate to zero, one obtains 
that the energy stored in the elastic solid is 



1 AfYf ,2 



/ eVLp CikCik + 



[1 + Ois')] 



2 Lp \ 1 — 2(Ts 

(A19) 

where V is the volume integrated over (assumed finite). 
In otherwords, any energy stored in the elastic solid that 
is due to the interaction between fibers, is smaller 
than the leading-order contribution which itself can be 
considered small. The leading order contribution only 
depends on the average number of broken fibers and thus 
does not alter the analytical form of Eq. ljA14|l . Thus, 
the energy stored in the elastic soHd plays no essential 
role in the damage model. 
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